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of Critical State Dynamics in Superconductive Films 



o 
o 

(N 



Yu. E. Kuzovle^Q 

Donetsk Physics and Technology Institute of NASU, 83114 Donetsk, Ukraine 

An one-dimensional model of magnetic flux penetration into thin strip-like superconductive fllm 
is subject to numerical analysis, which combines explicit account for specific oblate geometry of 
magnetic field lines around the film and a simplest phenomenology of the flux flow resistance under 
rigid pinning of vortices with temperature-dependent critical current. 
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1. During last decade, the classical theory of ther- 
momagnetic instability of excited critical states in type 
II superconductors [ij was supplemented with new nu- 
meric models and simulations 13 and new theory be- 
yond them 0, 0, . Both well explain the complicated 
dendritic fragmentation of magnetic flux observed exper- 
imentally 0- Especially impressive are simulations in 
terms of "atoms" with long-range interactions specific 
for vortices in films (se e also references and pictures 
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It seems interesting to compare so advanced "micro- 
scopic" approach with primitive "hydrodynamical" mod- 
els. Below we consider, perhaps, most simple one: an 
1-D model of magnetic flux penetration into strip with 
strong pinning of vortices. Our main issue will be explicit 
account for effects of real film's geometry, which involve 
non-uniformity of interaction between vortices, that is its 
long-range dependence on their absolute, in addition to 
relative, positions. Note that such non-uniformity may 
be inessential in respect to fine structure of captured 
magnetic flux (visualized by magneto-optics 0, Q) but 
essential in respect to its global characteristics (measured 
with SQUID or pickup coil HQ). 

2. A simple variant of the classical theory Jj reduces 
to equations (in CGS units) 
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J ^ — V X B 
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with magnetic inductance B , electric field E , electric 
current density J , temperature T , velocity of light c , 
specific heat C , "constitutive law" E = E{J, B, T) re- 
sulting from viscous dynamics of vortices and their pin- 
ning, and the dots deputizing for heat transfer. 

In case of superconductive film, one also needs in the 
Laplace equation for magnetic field out of film, in order 
to state a relation between tangential and normal compo- 
nents of the field at film's surface. At present, analytical 



solutions to this problem are known for two film shapes 
only: strip 0, and disk 0, ^| . The first case is sim- 
pler, therefore let our film will be a strip oriented along 
y-axis, with section \X\ < L/2, \Z\ < d/2, where L 
and d are film's width and thickness, respectively. 

Firstly, let us transform Eas lII2l into dimensionless 
form, choosing xq — L/2 to be unit of length (then 
the section becomes \X\ < 1, \Z\ < 6 = d/L), to — 
Xf)Hc2/ cpnJc as unit of time, with Hc2 being upper crit- 
ical field, pn - characteristic normal specific resistance, 
and Jc - maximum critical current density, jo = Jc as 
unit of current density, = 2'KdJclc as unit of mag- 
netic field, and eo = boxo/cto as unit of electric field. 

Secondly, carry through spatial averaging of B , J 
and other patterns along strip length, designating this 
operation by {. . .)y , and introduce the quantity 
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Applying (. . .)y to EqsQland using the Clem et al. 
and Brandt results .1^] it easy to obtain 
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du , (5) 
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where B — B{x,t) is normal component of {B)y addi- 
tionally averaged over thickness of the film, E ~ E{x, t) 
is similarly twice averaged K-component of E , and 
is external bias field (assumed uniform and perpendicular 
to film). The Eqs0| and |31 are applicable irrespective of 
any "dendritic" or other longitudinal non-uniformities, if 
only B possesses statistical uniformity. 

The integral in ((SJ is nothing but exact average tan- 
gential component of B at upper film's surface in the 
formal limit (5 — > . At that, we will keep in mind that 
physically our film is "not too thin" : d> \ with A being 
characteristic penetration depth. This allows to interpret 
magnetic flux density B{x,t) as direct measure of con- 
centration of vortices. 

What is for the averaging of the "constitutive law", 
E = E{J, B, T) , it is less trivial, since its result, some 
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effective law E = E{I,B,T) for averaged variables, can 
be sensible to any non-uniformities. At present, lacking 
something better, we will try something from traditional 
model laws. For instance (in our dimcnsionlcss units), 



E{I,B,T) = p{B)Ul,T) , 



In{I,T)= { 

I + Ic{T) 



p{B) = \B\ , h{T) = 1 - Q 



/ > Ic{T) , 
\I\<Ic{T) , 
I < ~Ic{T) , 



Tc-To 
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(8) 



Here p{B) represents the "flux flow resistance", Ic is 
critical value of J as a function of local temperature, 
and To and Tc are background (thermostat) and critical 
temperature, respectively. 

Thus the model assumes "rigid" pinning without creep, 
when vortices can move under |/| > Ic{T) only. At 
that, the excess current /„ = /„(/, T) directly deter- 
mines local drift velocity of vortices while the electric 
field E — \B\In local flux flow. Frequently J„ is treated 
conversely as "normal current" caused by motion of vor- 
tices. But we will see below that sometimes /„ turns 
into Meissner super-current shielding a part of film. In 
such "quasi-Meissner" states, in addition to drift of vor- 
tices, their diffusion may be important. We will take it 
into account if replace Eq.^] by 
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FIG. 1: Hysteresis loops of mean magnetic inductance, 
B{t) = J B{x,t) dx / L (black), and magnetization, M{t) = 
B{t) — Ho{t) (red), after start from ZFC (zero field cooled) 
state under saw- like Ho{t) with period 200, at parameters 
p = 10 , 7 = 2 and D ^ 0.03 . 



To finish writing of the temperature equation, include 
there heat diffusion across film and local heat exchange 
among it and its substrate. In the dimensionless form. 
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(13) 



A(/,T) 



, |/|</c(r) 

A , |/| > Ic{T) 



(10) 



so that diffusion occurs only when vortices get free from 
pinning. In view of the Einstein relation, the diffusivity 
A is definite function of the foregoing parameters. In our 
dimensionless units, A = 2T/fcL, with fc = ^oJcd/c 
being characteristic critical Lorentz force and pinning 
force acting per vortex. 

Boundary conditions to the diffusive term in Eq|51will 
state continuity of average normal component Bz of 
magnetic inductance at film's edges a; = ±1 . Outside 
of the strip it can be obtained from (see Q) 



B, 
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-1 



1 ~ n (u — w) 



where w = x + i\z\ , time is omitted, and Bz — 
{Bz{x,y, z))y ■ Particularly, in film's plane (at > 1 ) 



At very edges, of course, factor 1/y/x^ — 1 should be 
reasonably cut, e.g. by 1/V2d. 



p = 27TLdJ^/[c\T,-To)C] 



(14) 



(for simplicity, C and D assumed temperature inde- 
pendent). Estimates of dimensionless thermal diffusivity, 
thermal relaxation rate. Joule heating and vortex diffu- 
sivity for realistic films typically yield I?~0.01, 7~1, 
p 10, and A ~ 10"^ (ahhough values 10 100 times 
greater or smaller also may be realistic). 

3. The system of equations (I5|)- ((T^ was numeri- 
cally investigated using the third-order adaptive Runge- 
Kutta algorithm. The integral in Q was reduced to the 
Hilbert transform (at finite interval) in its turn performed 
through fast Fourier transform. It was found that sub- 
ject to values of p/j and \dHo{t)/dt\ the model has rich 
variety of regular and chaotic regimes of magnetic flux 
entry into film or departure from it. But we will not list 
all that. Our aim is to list only main properties of the 
system manifested in all regimes, and illustrate them by 
intermediate quasi-regular regimes with p/7 ~ 10 and 
dHoit)/dt^±0.1. 

Fig.l demonstrates clear signs of many magnetic flux 
avalanches with different magnitudes and durations. Rel- 
atively large-scale and reproducible avalanches alternate 
with quiet flux creep which on closer examination con- 
sists of many small-scale avalanches. Most tiny of them 
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FIG. 2: Two variants of total film's magnetization, M — 
'B ~ Ho (black) and M' introduced by Eas ll5ll6l Cred). via 
time, under the same conditions as in Fig.l. 

are indistinguishable from random computational errors 
which thus serve as mathematical noise generator imi- 
tating physical thermal fluctuations. Under the adaptive 
algorithm, a level of this noise occurs rather constant. 

From thermodynamical point of view, the mean flux 
density B~ J_j^^^B{x,t) dx/L and corresponding mag- 
netization M = B — Ho shown in Fig. 1 seem most nat- 
ural global magnetic characteristics of the film. But, ac- 
cording to from the point of view of any pickup coil 
the film is presented by a modified magnetization, 

nL/2 

M' = / W{x)[B{x,t) - Ho]dx , (15) 

J-L/2 

with W{x) being some weight function (normalized to 
unit). In particular, if a coil is placed far from film, then 

W{x) = -y/l-x^ (16) 

TT 

(in the dimensionless form). Fig. 2 shows time evolvent 
of such defined M' in comparison with M . 

It is popular to search for scale invariance in size dis- 
tributions of the avalanches. Then one firstly should in- 
troduce a criterion for their beginnings and endings. For 
example, let us define the borders between subsequent 
flux jumps as those inflection points of the plot B{t) 
where \dB(t)/dt\ has local minimums, and interpret in- 
crements or decrements of B{t) between these points 
(after multiplying them by film's width) as heights d^ 
of the jumps. Their probability density Wj{d^) (nor- 
malized histogram) is shown in Fig. 3. By request, one 
can find here good signs of the scaling. But dominant 
smooth part of this probability density corresponds to 
mini- avalanches invisible in Figs.l and 2. 
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FIG. 3: Histogram of magnetic flux jumps determined from 
"horizontal inflections" of B{t) (see body text), under the 
same conditions as in Figs. 1-2, in double logarithmic and (in- 
set) linear scales. 
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FIG. 4: Magnetic inductance at a; = 0.02 (black), a; = 0.1 
(red), X = 0.4 (blue) and x — 0.8 (green), as a function of 
external field rising linearly from ZFC state ( Ho{t) — O.OSt ). 
The gray straight line guides B = Ho , and x means distance 
from left film's edge. 

4. Consider peculiarities of the model which are defi- 
nitely caused by its flat thin-fllm geometry. 

First, even in simplest situation, when external fleld 
slowly and monotonously (linearly) goes up from ZFC 
state, local magnetic flux density inside fllm can be 
strongly non-monotone function of both time and po- 
sition. Second, in some vicinity of fllm's edges the flux 
density can significantly exceed external field. The latter 
circumstance as well as the temporal non-monotony are 
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FIG. 5: Typical instant spatial distributions of magnetic flux 
density B(x,t) (black), current I{x,t) (red) and critical cur- 
rent ±/c(x,t) (gray) across strip's width, at beginning stage 
of the process shown in Fig. 4. As in Fig. 4, x is dimensionless 
coordinate counted from left film's edge. 

demonstrated at Fig. 4 (in fact it relates to initial stage 
of the same process as in Figs. 1-3). 

Of course, the reason for such behavior is strong com- 
pression of magnetic field lines near film's edges (formally, 
at 6 —^ magnetic inductance at very edges can be ar- 
bitrary large). Sufficiently big avalanche drops the com- 
pression down, but then pinning enforces it to grow again. 

Example of spatial non- monotony is presented in Fig. 5. 
This snapshot is made just after avalanching. It is seen 
that the current became smaller then critical current 
and thus flux motion is stopped everywhere except the 
film's periphery, where slow creep (formed by local mini- 
avalanches) takes place. Note that in accordance with 
corresponding temperature pattern can be viewed if 
lift the lower gray curve —Ic{x) by unit. 

Clearly, such picture of B{x) as in Fig. 5 never could 
be observed (in similar process) under 3-D cylindric ge- 
ometry, since non- monotony of B{x) would mean sign 
reversal of J{x) (x —dB{x)/dx and thus reverse flux 
flow. But under film geometry, when differential relation 
between magnetic inductance and current changes into 
integral relation (O, the B(a;)'s non-monotony does not 
necessarily involve /(a;)'s sign reversal. Consequently, 
2-D geometry produces much more intricate patterns. 

For comparison. Fig. 6 shows what the model gives 
when thermomagnetic instability is out of play, under 
p/7 <C 1 . Then penetration of magnetic flux realizes 
merely as smooth drift, additionally stimulated by weak 
diffusion and noise when i?o(i) crosses zero (see below). 
At that, instant B{x,tys and /(x,t)'s patterns almost 
do not differ from ones at static critical state under the 
same but fixed external field. 

One more interesting thing seen from Fig. 4 is dis- 
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FIG. 6: Spatial distributions of magnetic flux density B{x, t) 
(black) and current I{x,t) (red) under linearly rising external 
fleld in absence of thermomagnetic instability (p/7 ^ 1 ). 
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FIG. 7: Local flux density in the middle part of the film at 
X — 0.5 (black) and in its center at 2; = 1 (red) under the 
same conditions as in Fig.l, with x being distance from left 
film's edge. 



cretization of local fiux density values deep in the film. 
Hysteresis loops drawn by them are shown in Fig. 7. 
These pictures prove that almost all avalanches begin at 
periphery of the film and therefore flux filling of its inte- 
rior comes from most powerful avalanches, while smaller 
ones never reach its middle. 

The evident presence of a large amount of randomness 
(partly chaos and partly noise) in evolution of B{x,t) 
and I(x, t) distributions inevitably results in breakdown 
of mirror symmetry of B{x,t) and anti-symmetry of 
I{x,t) , which is illustrated by Fig. 8. The breakdown oc- 
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FIG. 8: 

Typical picture of the mirror symmetry breakdown, pho- 
tographed at time 42.5 (when Ho = 3.4 ) during the same 
process as in Fig.l and Fig. 2. 
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FIG. 9: 

Top: Patterns of magnetic flux density B{x,t) (black) and 
current I{x,t) (red) during quasi-Meissner state before explo- 
sive flux reversal (thick curves) and quasi-critical state after 
it (thin curves). The gray lines show ±/c . 
Bottom: Distributions of flux flow Bin in the two above 
mentioned situations. For details see body text. 
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FIG. 10: Hysteresis loops of magnetization under the same 
parameters as in Figs. 1-2 but with 30 times greater (red curve) 
and 1000 times greater (black curve) A value. 

curs soon after start from ZFC and then never completely 
disappears, although the symmetry partly restores at 
borders between quadrants of the Hq — M diagram. 

5. Next, we focus our attention on two (mutually in- 
verse) biggest peaks of the magnetization in Figs.l and 
2 (at time moments t w 110 and t « 210) and also on 
red curve in Fig. 10. They represent not usual avalanches 
but flux reversal in the film when Ho(t) passes through 
zero. This process is illustrated also by Fig. 9. 

At that time (if saying about positive peak) earlier ac- 
cumulated positive magnetic flux decreases down to zero 
and then changes its sign. This realizes in two different 
ways: at first by flow of positive flux out of film, when 
vortices leave it, and afterwards by flow of negative flux 
into fllm, when anti- vortices enter it and annihilate with 
vortices in its interior. Obviously, transition from the 
"vortex departure stage" to the "anti- vortex entry stage" 
occurs when (i) magnetic field Bf''^^ at outward vicinity 
of film's edges becomes negative and then (ii) diffusive 
fiow of anti- vortices at film edges, Fav , becomes exceed- 
ing drift outflow of vortices, Fy = B{x = ±l)/„(a; = ±1) 
(importantly, anti-vortices do not contribute to drift be- 
cause their concentration is zero while the first "depar- 
ture stage" still lasts). 

According to EqE|(with the remark under it), 

Bfoe 1 r l^lBix)-Ho]dx (17) 

7rV2^ J-i V 1 T 2^ 

where upper sign corresponds to right-hand edge. This 
expression clearly implies that B^'^^'^ crosses zero before 
HQ{t) does it. The diffusive anti-flux flow is proportional 
to gradient of Z-component of magnetic inductance at 
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the very edges and can be written as 

F,, « A [Bix = ±1) - Bf9-]/S 



(18) 



by definitions of the appearing quantities (aU smoothed 
over film thickness). The EasllTI and 1181 lead to estimate 



A 



1 ± x „ , , dx 

B{x)—= 

iTx ttV2 



(19) 



valid when Ho{t) already crossed zero, at Hq < (of 
course, if the transition not yet happened, and B{x) 
looks as in Fig. 9). 

Hence, at sufficiently small A/ 6^^'^ the transition from 
"departure stage" to "entry stage" can be strongly de- 
layed. Just such example is presented by Figs. 1-2, where 
the annihilation switches on at t w 110 only, when 
Hq « —0.8 (i.e. under Hq comparable with full pen- 
etration field!). 

This results in strong "supersaturation" of the vortex 
system and, as consequence, fast explosive transition, 
which manifests itself in Figs. 1-2 as high sharp steep in 
B{t) and M{t) plots. Increase of A , under the same 
thickness, shortens "departure stage" but prolongs and 
softens "entry stage" . Correspondingly, the steep in mag- 
netization curve decreases and at sufficiently large A it 
disappears at all, as Fig. 10 shows. 

Of special interest is spatial structure of the depar- 
ture stage. As top of Fig. 9 shows, due to the outflow 
of vortices their peripheral concentration drops to very 
small values. However, the outflow continues. Therefore, 
E = \B\In must be approximately independent on spa- 
tial coordinate X . Indeed, thick blue curve at bottom 
of Fig. 9 conflrms good satisfaction of this requirement, 
while thick red curve at top of Fig. 9 shows that it is 
satisfied owing to strong growth of "normal current" /„ 
close to film's edges. There ~ Ic or even |/„| ::§> Ic , 



and this current well screens film's periphery. Hence, we 
see something like Meissner super-current and Meissner 
state! 

In essence, of course, this "quasi-Meissner" state rep- 
resents not a quasi-static state as true Meissner states, 
but sooner quasi-stationary dissipative structure, just be- 
cause it is maintained by motion of vortices (up to a mo- 
ment when it will be destroyed by birth of anti- vortices). 
In this sense, /„ by right can be qualified as dissipative 
"normal" current. 

Thin curves at Fig. 9 present distributions of B{x, t) 
and I{x, t) just after the explosive annihilation and flux 
reversal. Evidently, annihilation still not finished at 
film's center, but in the rest of film |/| w /c ~ 1 , i.e. 
negative flux evenly fills it. Corresponding plot of flux 
flow E = \B\In is shown by thin curve at bottom of 
Fig. 9 (we plot Bin instead of \B\In in order to un- 
derline transition from vortex flow to anti- vortex one). 
Thus, we observe typical "quasi-critical" state (but next 
avalanches are not far off). 

6. To resume, we tested a simple model of magnetic 
flux penetration into films of type II superconductor with 
temperature-sensible pinning. Advantage of the model is 
that it is based on exact extremely non-local relations be- 
tween fiux density distribution and current distribution 
in film. Defects of the model are in rather rough phe- 
nomenology of dissipative and thermal processes. The 
constituents of the model are not original, but, to the 
best of my knowledge, the model as a whole still was not 
under careful investigation. I tried to demonstrate that 
it is enough substantial and interesting and even can pro- 
duce useful tips for understanding experimental data. 

I am very grateful to Dr. Yu. Medvedev and Dr. 
V. Khokhlov for fruitful information and attracting my 
attention to the subject of this paper. 
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